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In this paper we perform a detailed analysis of the effect of various approximations which have 
been used in the literature to compute the detection rates of compact binary coalescences for in- 
terferometric gravitational wave detectors. We evaluate the detection rates for the coalescence of 
BH-BH, NS-NS, and BH-NS binaries taking into account their specific statistical properties ob- 
tained from population synthesis models, the cosmic star formation rate history, and the effects of 
redshift on the emitted gravitational wave signals. The results are compared with those obtained 
with procedures that are based on different levels of approximations, such as adopting averaged 
values for the total mass and symmetric mass ratio for all the systems of a binary population, using 
these to compute the horizon distance for individual detectors, or estimating the coalescence rate 
density within this distance by its local value. We show that these approximations introduce a bias 
which depends not only on the approximation which is used, but also on the statistical properties 
of the considered binary population. 
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I. INTRODUCTION 



The second generation of gravitational wave interfer- 
ometers, Advanced Virgo^ and LIGO^, will push ten 
times farther their observational reach, opening a much 
wider region of the Universe to observations. Fur- 
thermore, the third-generation ground-based observatory 
Einstein Telescope (ET)^, which is currently in a design 
study phase, is expected to improve the sensitivity by 
a further factor of ten [T^ T . This will allow to detect 
gravitational wave (GW) signals emitted by astrophysi- 
cal sources, for instance the coalescence of compact bi- 
naries in which we are primarily interested, with suffi- 
ciently high signal-to-noise ratios (SNR) and rates, to 
allow precise astronomical studies of these sources [3] . In 
addition, the high sensitivity of these detectors will al- 
low to see beyond our local universe, accessing regions of 
redshift where the role of the star formation history and 
of cosmological models becomes important, and could be 
tested ([IHB] and references therein). 

The detection of a GW signal is based on the corre- 
lation of the output data stream of the detector with 
an expected signal, which is typically in the form of a 
template. Given the high number of parameters needed 
to model a source, a complete survey of the parameter 
space is often impossible. For this reason it is important 
to know what are the systems that are more likely to be 
detected and their expected properties [II [8] . In recent 



years several papers have been published where the coa- 
lescence rates of compact binaries, and the corresponding 
detection rates by GW detectors, have been estimated by 
using several approximations ^9r.22j . 

In this paper we study in detail the effects of these 
approximations on the detection rates of advanced de- 
tectors, also showing to what extent these effects de- 
pend on the binary population model which is used. We 
choose, as an example, two models generated by the code 
St art rack, and publicly available^. Details on the phys- 
ical inputs of the simulation are given in [23 . The binary 
populations are the result of the evolution of a large sam- 
ple {N = 2 • 10^) of "zero- age" binaries. 

Using these data, we compute the coalescence and de- 
tection rates of binaries composed of two black holes, 
BH-BH, two neutron stars, NS-NS, a black hole and a 
neutron star, BH-NS, using a completely consistent pro- 
cedure based on the appropriate use of the statistical 
properties of the populations, and on the redshift depen- 
dence of the star formation rate history [TOl |Tl] |24l |25] . 
The detection rates are given for Advanced Virgo and 
LIGO, and for ET. 

These rates are then compared with those computed 
by applying, to the same populations, the different ap- 
proximations used in the literature, to evaluate the bias 
they introduce (if any), to show how this varies for the 
Advanced detectors in different configurations, and to 
understand its dependence on the population properties. 

The statistical properties of coalescing binary popu- 
lations depend on the physical inputs of the population 
synthesis codes, which may vary significantly from one 
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model to another. The choice of the two models we con- 
sider in this paper, together with their main statistical 
properties, are briefly described in section [TH 

We will show that in order to obtain detection rate 
estimates which depend only on the properties of the 
binary population, we need to use a fully consistent ap- 
proach, since each of the approximations used in the liter- 
ature leads to over/underestimate these rates by amounts 
which depend on the population properties. 

The plan of the paper is the following. In section [ll| we 
describe the fully consistent procedure to compute the co- 
alescence and detection rates for a simulated population 
of compact binaries. In section |III| we describe various 



approximations used to evaluate the rates [H [20l |2T] , ap- 
ply them to the two population models we consider, and 
compare the results with those found with no approxi- 



mation. Concluding remarks are given in section IV 



II. RATE COMPUTATION PROCEDURE 

The binary population models we use to compute co- 
alescence and detection rates, are chosen from the Stan- 
dard model of [23] , since this model is the flducial one, be- 
ing based on conservative assumptions. We consider the 
two submodels A which correspond to progenitor stars 
with solar and subsolar metallicity. The corresponding 
populations exhibit different mass distributions, and this 
is the feature which is of major interest for us. Here- 
after we will refer to these two models as model A (solar 
metallicity) and model B (subsolar metallicity). 

Here we report some statistical properties of each 
model, which will be useful in the following to under- 
stand the results. First of all, it should be mentioned 
that, as a result of different physical inputs (which we 
do not discuss here), the number of binaries coalescing 
within the Hubble time, differs in the two models: for 
model A there are 1753 BH-BH binaries, 5084 NS-NS and 
345 BH-NS; for model B the population is made of 15650 
BH-BH, 1767 NS-NS and 769 BH-NS binaries. Thus, the 
lower metallicity used in model B causes an increase of 
the number of BHs that are formed and of their mean 
mass as well. 

The average total mass M and symmetric mass ratio 
fj are: 
for model A 



(1) 



M = 


= 15.41Mo f? = 


= 0.248 for BH-BH 


M -- 


= 2.43M0 f? = 


= 0.248 for NS-NS 


M -- 


= 9.94Mq f? = 


= 0.150 for BH-NS 


model B 







M = 3O.66M0 f] = 0.246 for BH-BH (2) 

M = 2.5OM0 f] = 0.249 for NS-NS 
M = 11.72M0 f] = 0.117 for BH-NS 

The procedure we adopt to compute the detection rates 
from the population synthesis statistics is the following. 



Starting from the star formation rate density at a given 
redshift, Pi,{z) (in units of Mq yr~^Mpc~^), we derive 
the binary progenitors birth rate per comoving volume 
(in units of yr~^ Mpc~^) as. 



where 



dR 
dtdV 



A: 



(z) = Ax p^{z) 



/bin ^ /si 



2{m*) ' 



(3) 



(4) 



/bin is the binarity fraction, (m^) is the average stellar 
mass and /sim is the fraction of binaries simulated by the 
population synthesis code. The latter quantity accounts 
for the fact that the range of progenitor masses used in 
the simulation is different from the range of validity of 
the adopted IMF. The factors /sim and (m^) depend on 
the adopted stellar Initial Mass Function (IMF), <l>(m), 
as 



{m^) 



Jg dm^{m) 



Jq ^ dmm^{m) 



(5) 



(6) 



In [2T the authors set /bin = 2/3 and initialize the stellar 
progenitors in the mass range [5 — 15O]M0 according to 
a Kroupa IMF. 

We extract from the simulated sample of Ntot binaries, 
those systems with total mass G (M, M-\-dM)^ symmetric 
mass ratio G {r]^r] -\- drf)^ and delay time G {td^ td + dtd), 
which form BH-BH, NS-NS, and BH-NS systems; the 
delay time is deflned as the time interval between the 
formation of the binary progenitor system and its coa- 
lescence. Their number is Ni{M^r]^td)^ where i refers to 
BH-BH, NS-NS, BH-NS systems. 

The number of binaries coalescing at redshift Zc, with 
total mass M and symmetric mass ratio 77 , per unit of 
M, T] and comoving volume is (Mpc~^yr~^MQ^?7~^): 



rcoa>(M,,,.e) = / ^-P^(^f)^^(^;^^^^) du ; (7) 



l + Zf 



N^ 



tot 



the redshifts of progenitor binary formation, Zf, and of 
binary coalescence, Zc, are related by the delay time as 
follows 



td 



«/ 



dz , 



and 



Ho{l + z)^^m{l + zf + {l-nm) ■ 



(8) 



(9) 



In this work we assume Hq = 73 km~ Mpc~^ and 
Qrn = 0.24, consistently with the choice done in 06] to 
derive our adopted cosmic star formation history, which 
reproduces the observational data available at 2; < 8 [27]. 
The coalescence rate density pc{zc) (Mpc~ yr~^), i.e. 
the number of compact binaries coalescing at redshift Zc, 
per unit time and unit comoving volume, is 



and for model B 



Pc{zc) = rcoai{M,r],Zc)dMdr] 



(10) 
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FIG. 1: (color online) The coalescence rate density computed 
with Eq. (10), is plotted versus redshift, respectively for BH- 
BH, NS-NS and BH-NS binaries, for model A and B. 

A (upper panel) and model B (lower panel). We define 
the value p|°^^^ = pc{zc = 0) as the "local merger rate 
per unit comoving volume" . For model A we find 

^locai ^ 4 54 X iQ-2 Mpc-^ X Myr-^ (BH-BH)(11) 
^locai ^ ;l 43 X iQ-i Mpc-^ X Myr-^ (NS-NS) 



•local 



1.40 X 10-2 Mpc"^ X Myr"^ (BH-NS). 



•local 
•local 



4.32 X 10-^ Mpc"^ X Myr"^ (BH-BH)(12) 
4.69 X 10-2 Mpc"^ X Myr-^ (NS-NS) 
^locai ^ 2.73 X 10-2 Mpc-^ X Myr-^ (BH-NS). 

To compute the detection rate, we sum over the sources 
which are within a given interferometer reach. Assum- 
ing that the signals are extracted from noise using the 
standard matched filtering technique, the Signal to Noise 
Ratio (SNR), averaged over sky positions and angles is 



{SNR) 



2 3 r2 (l + z)2 roo 



df- 



1 dE, 



2 {5^^dL{z)^Jo ' pShU) df, 



[(1 + ^)/]] 

(13) 

where / is the observed frequency, related to the source 
emission frequency /e by 



/ = 



/e 



l + Z 



(14) 



dEe/dfe is the source energy spectrum, and Sh is the 
detector noise spectral density. For BH-BH systems we 
shall use the most recent hybrid waveforms produced by 
numerical relativity groups [28 . These describe the full 
inspiral-merger-ringdown phases of coalescence, account- 
ing also for non-precessing spins. However, since the 
simulated populations we use do not include spins, we 
shall assume that all BH-BH systems have zero spin. For 
NS-NS and BH-NS binaries we use the chirp waveform, 
truncated to the ISCO. Moreover we will consider vari- 
ous noise curves of advanced detectors, namely Advanced 
LIGO (ALIGO) in four of its planned configurations 
[29]^, namely zero detuning high laser power ( ALIGO 1), 
optimized for 30-30 solar mass bh binary (ALIG02), no 
signal recycling mirror (ALIG03), and optimized for NS- 
NS binaries (ALIG04), Advanced Virgo^ ( AVirgo) , and 
ET'^ in three different configurations namely ET_B, ET_C 
and ET_D. 

We define Zmax{M^r]) as the maximum redshift at 
which a system with assigned (M, 77) is detectable 
with SNR=8, averaging over position and orientation; 
ZmaxiM^ 77) is obtained by inverting Eq. (fl3]). D^g,^{M^ 77) 
is the corresponding luminosity distance. 

Moreover, we define Zh(M, 7^) and D\^{M^r])^ as the 
horizon redshift and horizon distance, computed as be- 
fore, but assuming optimal orientation. 

The detection rate per unit total mass M and mass ra- 
tio 7^, Nm.ti'! expressed in units [yr"^M0^77"^], is found by 
integrating the function r'coai(^5 Vi ^c) up to z^q,^{M^ 77), 
as follows: 



5 see LIGO document Iigo-t0900288iv3 

6 wwwcascina.virgo.infn.it/ advirgo/docs/ AdV_refsens_100512.txt 
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No Approximation detection rates A^det [yr ^] 


detector 


BH-BH 


NS-NS 


BH-NS 


model A 


model B 


model A 


model B 


model A 


model B 


AVirgo 


73 


5.67 • 10^ 


2.3 


0.82 


3.7 


6.61 


ALIGOl 


203 


1.52-10^ 


6.0 


2.10 


9.2 


17.8 


ALIG02 


134 


1.07-10^ 


3.9 


1.35 


6.2 


12.0 


ALIG03 


105 


8.34 • 10^ 


3.1 


1.08 


5.0 


9.60 


ALIG04 


254 


1.84-10^ 


7.6 


2.66 


11.8 


22.0 


ET_B 


1.80-10^ 


1.56-10^ 


2.65 • 10^ 


8.01 • 10^ 


1.79 • 10^ 


3.62 • 10^ 


ET_C 


1.84-10^ 


1.59-10^ 


4.48 • lO'' 


1.35 • lO'' 


2.47 • lO'' 


5.39 • lO'' 


ET_D 


1.81 • 10^ 


1.57-10^ 


2.64 • 10^ 


7.95 • 10^ 


1.87-10^ 


3.86 • 10^ 



TABLE L Detection rates computed from Eq. (17), for the simulated populations of BH-BH, NS-NS and BH-NS (model A and 
B) and various detector sensitivities (first column) (see text). 



r2;max(M,77) ^Y 

NM,r^{M,r])= I f^oAM,r],Zc)j^dZc, (15) 



where the function 
dV 



dz 1 -\- z 



(16) 



is the comoving volume element and dL{z) is the lumi- 
nosity distance. It has to be noted that Nm,?] depends 
on M and 77 , not only through the coalescence rate rcoab 
but also through the integration limit. This means that, 
for each system (i.e. for each value of M and r] ), the 
coalescence rate is integrated up to a different distance. 
The integrated detection rate, i.e. the final output of 
our procedure, in units [yr~^] is: 



iVdet = 



f dM f NM,^{M,r])dr] 

f f r^— (- 

dM df] 



(17) 

rco^\{M,r],Zc)-^dZc . 
dzc 



In the following we will refer to the results obtained with 
this procedure as the "no-approx" rates. 

In Table [l| we show the detection rates obtained for 
the three considered families of compact binaries, using 
model A and B. The uncertainty in these results, due to 
numerical integration, is of a few percent. 

The data show that, for model A and B, the highest 
rates are attained by BH-BH systems, both for advanced 
and for third generation detectors. 

As already noted in [20], where detection rates where 
computed only for advanced detectors, including metal- 
licity effects in stellar evolution results in a much higher 
detection chance for BH-BH and BH-NS coalescing bi- 
naries. This is because in model B the number of these 
systems is larger than in model A, and the average mass 
of systems including BHs is higher (see Eqns. (fTl) and ([2| 
and the text above). 



Horizon distances Dh [Mpc] 


detector 


model A 


model B 


BH-BH 


NS-NS 


BH-NS 


BH-BH 


NS-NS 


BH-NS 


AVirgo 
ALIGOl 
ALIG02 
ALIG03 
ALIG04 


1980 
2980 
2520 
2290 
3260 


371 
518 
445 
412 
562 


959 
1390 
1200 
1100 
1510 


4030 
6110 
5210 
4580 
6770 


381 
533 

458 
423 

578 


955 
1390 
1200 
1100 
1500 



TABLE II: Horizon distances of advanced detectors (in Mpc), 
computed for the mean values of mass and mass ratio given 
in Eqns. (IT]) and ([2I, for the three families of binaries and for 
models A and B. 



Finally, in Table [llj we give the values of the horizon 
distance, D\^ = D\^{M^ ff)^ which corresponds to the mean 
values of mass and mass ratio given in Eqns. ([l]) and ([2|; 
this quantity will be useful in the following discussion. 



III. COMPARISON WITH OTHER 
PROCEDURES 



Various approaches have been used in the literature to 
compute detection rates, starting from the coalescence 
rates produced by population synthesis simulations [9l 
[T2H22] . A common feature of these approaches is that 
the detection rate is evaluated by multiplying the merger 
rate by a properly defined detection volume 14, rather 
than using Eq. (17). 



In [9] the authors adopt the following procedure. As- 
suming that all black holes have M = 10 Mq and all 
neutron stars have M = 1.4 M©, they define the aver- 
age horizon distance D\^^ such that a binary system is 
detected with a SNR > 8 by a given detector. 

A detection volume is defined as the Euclidean volume 



Approximation 1: detection rates [yr ^] 


detector 


model A 


model B 


BH-BH 


NS-NS 


BH-NS 


BH-BH 


NS-NS 


BH-NS 


NLt 


iVdet 


NLt 


iVdet 


NLt 


iVdet 


NLt 


iVdet 


nL 


A^det 


nL 


A^det 


AVirgo 


79 


73 


2.4 


2.3 


3.8 


3.7 


3.98 • 10^ 


5.67- 10^ 


0.78 


0.82 


7.86 


6.61 


ALIGOl 


210 


203 


6.2 


6.0 


9.3 


9.2 


1.07-10^ 


1.52 • 10^ 


2.01 


2.10 


21.0 


17.8 


ALIG02 


140 


134 


4.0 


3.9 


6.5 


6.2 


7.17-10^ 


1.07-10^ 


1.30 


1.35 


14.1 


12.0 


ALIG03 


107 


105 


3.2 


3.1 


5.3 


5.0 


5.36 • 10^ 


8.34 • 10^ 


1.04 


1.08 


11.3 


9.60 


ALIG04 


260 


254 


7.9 


7.6 


12.3 


11.8 


1.36 • 10^ 


1.84-10^ 


2.55 


2.66 


25.2 


22.0 



TABLE HI: The detection rate A^det? evaluated from Eq. (21 ) (Approx 1) for model A and B. A^det fo^ BH-BH systems (column 
2), is compared to the no-approx rates, A^det, computed from Eq. (17) (column 3). The additional columns refer to the same 
quantities evaluated for NS-NS and BH-NS binaries. 



within Dh- 



H 



3 \MpcJ 



(18) 



The detection rate is then computed by multiplying the 
Galactic coalescence rate, pf^^ expressed in units Myr~^ 
MWEG"^ (Milky Way Equivalent Galaxy), by the num- 
ber of MWEGs inside 14, G{Dy,), 



where 



iVdet = Pf ' X G{Dh), 



G(Dh) = VhX (0.0116) X 



2.263 ' 



(19) 



(20) 



The numerical factor 0.0116 is an estimate of the local 
number density of MWEGs in units Mpc~ obtained in 
[30 ; the factor 1/2.26^ is needed to average over sky 
position and orientation. This factor and the distance 
i^h in Eq. (18), in ref. [9 are computed neglecting the 
effect of redshift (see their footnote 90). 

The Galactic coalescence rate, p^^^ can be estimated 
following two different procedures, i.e. by using a con- 
stant Galactic star formation rate [3T1[32], or by normal- 
izing the number of binaries produced by a population 
synthesis code to the observed Galactic supernova rate 
[Tll|2T]. Using the former procedure, in [23 they find 



• P, 



,Gal 



^.2 MWEG"^ X Myr"^ for BH-BH, 



• p^^^ = 23.5 MWEG"^ X Myr"^ for NS-NS, 

• p^^^ = 1.6 MWEG"^ X Myr-^ for BH-NS. 
for model A and 

• p^^^ = 73.3 MWEG"^ X Myr"^ for BH-BH, 



.^Gal 



^.1 MWEG"^ X Myr-^ for NS-NS, 



p^^^ = 3.4 MWEG"^ X Myr"^ for BH-NS. 



for model B. If the detection rate is estimated using the 

local 

Mpc" 



coalescence rate 

"^ X Myr"^ , Eq. (Il9|) reads 



3er unit volume, p^ ^ in units 



iVde 



ilocal 



^lOCal y^Y^y^ 



2.263* 



(21) 



It should be noted that the ratio between the local and 
the Galactic rates depends on the adopted cosmic star 
formation history model, on the Galactic star formation 
rate, as well as on the statistical properties of the binary 
population (delay time distribution). In [9 a constant 
value of pI^^^YpI^^ = 0.0116, has been adopted for all 
compact binaries considered in their study. In summary, 
the approximations done in [9 to estimate the detection 
rates predicted by different binary populations and dif- 
ferent detectors are: 

1. Approx 1: use of a single horizon distance I^h for all 
systems of a given binary population, correspond- 
ing to the mean values of mass and mass ratio. 

2. Approx 2: evaluation of the detection rate by mul- 
tiplying the local value of the coalescence rate den- 
sity, pc{z = 0), by the Euclidean volume, within 



the horizon distance I^h, as in Eqns. (19) or (21). 



3. Approx 3: evaluation of the horizon distance ne- 
glecting the redshift effects included in Eqns. (13) 
and (Il4|). 



In what follows, we investigate the effects of these ap- 
proximations; since in [9j they are used all together, we 
will show the consequences of applying approximation 1 , 
1+2 and 1+2+3 . 



A. Approximation 1 

Approximation 1 consists in using the same horizon 
distance for all systems of a binary population, regardless 
of their individual total mass and symmetric mass ratio. 
To evaluate the effect of this approximation, we use the 
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FIG. 2: The total mass distribution for BH-BH systems for model A (left) and B (right) is plotted in the upper panels. In the 
lower panels we plot the integrand of Eq. (22) as a function of M and redshift Zc, using a gray scale proportional to its value. 
The dashed line corresponds to 2max(M, ry), the solid line corresponds to the redshift up to which the system with a given total 
mass and 77 = fy would be detected with an average SNR of 8. 



values of 5h given in Table [llj which correspond to the 
average values of M and rj given in Eqns. ^ and ([2]). 
We then compute the detection rate N^^^^ by integrating 
Eq. (17) up to the redshift Zmax 



nL 



where Zr. 



dM 



dZr 



dV 

rcoai{M,r],Zc)-r- df] , (22) 
azc 



ble 



III 



= Zmax{M ^fj). The results are shown in Ta- 
For model A, the difference between N^^^ and the 
no-approx rates A^det is of a few percent, for all binary 
populations, and N^^^ is always larger than A/'^et- 

For model B, conversely, the difference between N^^^ 
and A^det ranges between -26% (ALIG04) and -36% 
(ALIG03) for BH-BH, and between +13% (ALIG04) 
and +19% (AVirgo) for BH-NS. For NS-NS binaries the 
difference between N^^^ and A'det is, as for model A, of a 
few percent for all interferometers. Moreover, for BH-BH 
and NS-NS families, A/'^^t i^ systematically smaller than 
A'det, while it is always larger for BH-NS. 

These data show that the effect of Approx 1 depends 
on the model and on the kind of binary population. To 



clarify this dependence, which will affect also the results 
obtained with Approx 2 and 3, in the upper panel of 
Fig. [2] we plot, as an example, the total mass distribution 
for BH-BH systems for model A (left) and B (right). In 
the lower panels we plot the function 



rcoa\[M,r],Zc)-r- drj , 
aZc 



(23) 



which appear in Eq. ( [22| ), as a function of M and of the 
redshift Zc^ using a gray scale proportional to its value. 
In these panels, the horizontal dashed line corresponds to 
^max(^, ^), whereas the solid line represent the function 
^max(^7 ^)^ i-6- ^he redshift up to which the system with 
total mass M and symmetric mass ratio fj would be de- 
tected with an average SNR of 8. When the function (23) 
is integrated up to z^ to compute the rate as in Eq. (22), 
the contribution of masses lower than the crossing point 
between the two lines is overestimated, whereas that of 
masses above that point is underestimated, with respect 
to their contribution in Eq. (17). 

From Figure [2| we see that for model A, the two con- 
tributions nearly balance, whereas for model B, which 



Approx 1+2: detection rates [yr ^] 
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model B 


BH-BH 
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BH-BH 


NS-NS 


BH-NS 




iVdet 




iVdet 


^^det 


iVdet 
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^^det 


iVdet 


^^det 


^det 


AVirgo 


91 


73 


2.7 


2.3 


4.6 


3.7 


5.62 • 10^ 


5.67- 10^ 


0.94 


0.82 


9.20 


6.61 


ALIGOl 


273 


203 


6.9 


6.0 


12.1 


9.2 


1.73 • 10^ 


1.52 • 10^ 


2.46 


2.10 


24.7 


17.8 


ALIG02 


173 


134 


4.4 


3.9 


7.9 


6.2 


1.09 • 10^ 


1.07-10^ 


1.56 


1.35 


16.4 


12.0 


ALIG03 


134 


105 


3.5 


3.1 


6.3 


5.0 


8.09 • 10^ 


8.34 • 10^ 


1.25 


1.08 


13.0 


9.60 


ALIG04 


347 


254 


8.8 


7.6 


15.7 


11.8 


2.30 • 10^ 


1.84-10^ 


3.12 


2.66 


31.5 


22.0 



TABLE IV: The detection rate A^det, evaluated from Eq. (21 ) using Approx 1+2, for model A and B. A^^et for BH-BH systems 



(column 2), is compared to the no-approx rates, A^det, computed from Eq. (17) (column 3). The additional columns refer to 
the same quantities evaluated for NS-NS and BH-NS binaries. 



Horizon distances [Mpc] 


detector 


model A 


model B 


BH-BH 


NS-NS 


BH-NS 


BH-BH 


NS-NS 


BH-NS 


m 


5h 


m 


5h 


Dl 


5h 


Dl 


5h 


m 


5h 


Dl 


5h 


AVirgo 
ALIGOl 
ALIG02 
ALIG03 
ALIG04 


1540 
2120 
1860 
1720 
2280 


1980 
2980 
2520 
2290 
3260 


347 
474 
412 
383 
511 


371 
518 
445 
412 
562 


830 
1140 
1000 
931 
1230 


959 
1390 
1200 
1100 
1510 


2650 
3650 
3200 
2960 
3920 


4030 
6110 
5210 
4580 
6770 


356 
486 
422 
392 
524 


381 
533 

458 
423 

578 


829 
1140 
1000 
934 
1230 


955 
1390 
1200 
1100 
1500 



TABLE V: Average horizon distances (in Mpc) for the three families of binaries and for models A and B. Df^ is computed 
neglecting redshift contributions, whereas these are included in D^. Both quantities are evaluated using the mean values M 
and fy given in Eqns. (IT]) and (l2|. 



has a bimodal mass distribution, the high mass contri- 
bution, which is relevant, is neglected. This explains why, 
for the BH-BH population of model B, N^^^ is systemati- 
cally lower than A/'det- A similar explanation holds for the 
BH-NS population of model B, for which the low mass 
contributions is larger than the high mass one. Thus, 
the rate estimated using Approx 1 can either be larger or 
smaller than A^det by a different amount, depending on 
shape of the mass distribution. It has to be noted that 
the effect of this approximation is not related to the de- 
tector sensitivity because, for example, for model B the 
largest effect is for ALIO 03, which does not have the 
largest horizon distance (see Table [Tl|) . 



B. Approximation 1+2 

Approximation 2 assumes that the coalescence rate pc 
is constant through the whole detection volume Vh^ and 
equal to the value at z = (Eqns. ([l9| or (21)). As a 
consequence, the dependence of the coalescence rate on 
redshift is completely neglected. 

In addition, the detection rate is obtained by multi- 
plyin g pjf *^^^ = f)c{z = 0) by the Euclidean volume Vh 
(Eq. lis )) within I)^ax(M, 7/). 



To assess the bias introduced using these approxima- 
tions, we compute the detection rates, A/'^^f ^^• 



^det =3^ ^max 



(M,7))3xp|?^^i 



(24) 



The results are given in Table IV, where these are com- 
pared to the detection rates A/det computed in section [ll| 
with no approximation. 

Let us consider model A first. The values of A^^^^ 
are always larger than those of A'det, i-e. Approx 1+2 
overestimates the rates. The most significant differ- 
ences are for BH-BH, and vary from +24% (AVirgo) 
to +37% (ALIG04); for NS-NS they vary from +13% 
(ALIG02,ALIG03) to +17% (AVirgo) and for BH-NS 
from +24% (AVirgo) to +33% (ALIG04). Since, as dis- 
cussed above, for model A Approx 1 does not introduce 
a significant bias in the estimated rates, it is clear that 
the differences are due mainly to Approx 2 alone. 

These can be traced back to the contribution of the two 
assumptions: the first is that pc is assumed to be constant 
and equal to the z = value, while it is clear from Fig.[l] 
that this is not the case. This assumption alone would 
lead to an underestimate of the rates for all families be- 
cause, in the range of interest for z, the function pc{z) is 
always larger than the local value. Conversely, the second 



Approx 1+2+3: detection rates [yr ^] 


detector 


model A 


model B 


BH-BH 


NS-NS 


BH-NS 


BH-BH 


NS-NS 


BH-NS 


at! 23 


iVdet 




iVdet 


Arl23 
^^'det 


iVdet 


Arl23 
^^'det 


iVdet 


Arl23 
^^'det 


iVdet 


Arl23 
^^'det 


iVdet 


AVirgo 


63.3 


73 


2.4 


2.3 


3.7 


3.7 


3.0 • 10^ 


5.67 • 10^ 


0.85 


0.82 


7.43 


6.61 


ALIGOl 


166 


203 


6.1 


6.0 


9.0 


9.2 


7.85 • 10^ 


1.52-10^ 


2.15 


2.10 


18.4 


17.8 


ALIG02 


111 


134 


3.9 


3.9 


6.1 


6.2 


5.31 • 10^ 


1.07-10^ 


1.40 


1.35 


12.5 


12.0 


ALIG03 


88.9 


105 


3.2 


3.1 


4.9 


5.0 


4.21 • 10^ 


8.34 • 10^ 


1.12 


1.08 


10.1 


9.60 


ALIG04 


204 


254 


7.7 


7.6 


11.5 


11.8 


9.76 • 10^ 


1.84-10^ 


2.71 


2.66 


23.2 


22.0 



TABLE VI: The detection rates computed using Approx 1+2+3, i.e. neglecting the effect of redshift, and assuming that the 
merger rate is constant through the detection volume, and equal to its local value. The detection rates A^det, computed with 
the no-approx procedure (Eq. (17)) are given for comparison. 



assumption, i.e. that the local coalescence rate is multi- 
plied by the Euclidean volume, leads to an overestimate 
of the rates, because the Euclidean volume is larger than 
the cosmological one, and this effect is larger than the 
former, so that the combined result is an overestimate of 
the detection rates for all the three families. 

For model B, the detection rates N^^^ given in Ta- 
ble IV , differ (in percentual) from the corresponding val- 



ues of N^et as follows: for BH-BH from -3% (ALIG03) 
to +25% (ALIG04), for NS-NS from +15% (AVirgo) 
to +17% (ALIG01,ALIG04), for BH-NS from +35% 
(ALIG03) to +43% (ALIG04). By comparing N^^^ in 
table IIHI with N^^^ in table |IV| we see that also in this case 



the effect of Approx 2 alone is to raise the rates obtained 
with Approx 1 alone. It should be noted that, contrary 
to the results of Approx 1, the effect of Approx 2 depends 
on the reach of the interferometer, and is larger for in- 
terferometers with higher sensitivity. This is particularly 
evident for model A, where the effects of Approx 1 are 
negligible, while those of Approx 1+2 increase going from 
AVirgo to ALIG04, which has the highest sensitivity. 

For model B the effect of Approx 1 is not negligible, 
so that the effect of Approx 1+2 varies between interfer- 
ometers and binary families. 



C. Approximation 1+2+3 

Approximation 3 does not include redshift contribu- 
tions in the SNR evaluation. This will reflect in different 
values for Dh and I^max which we will denote as Df^ and 
^max- In Table 



ileMi 



we compare the horizon distances Dh 
given in Table |II[ and which include redshift contribution, 
Dl{M,fj). 

Horizon distances computed neglecting the redshift are 
smaller than those which include it, and this reflects into 
smaller detection rates, as can be seen by comparing the 
rates N^'^^ in table IV for approx 1+2, with N^^^, given 
in Table VI] and computed for Approx 1+2+3. As ex- 
pected, the difference is larger for interferometers with 
larger horizon distances. 



Turning to the comparison between N^^^ and A^det, 
shown in Table |VI[ we see that the percentage differ- 
ences for the three families of binaries range within: 
for model A, from -13% (AVirgo) to -18% (ALIGOl) 
for BH-BH, from (ALIG02) to +4% (AVirgo) for 
NS-NS, and from (AVirgo) to -3% (ALIG04) for 
BH-NS; for model B, from -47% (AVirgo) to -51% 
(ALIG02) for BH-BH, from +2% (ALIG01,ALIG04) to 
+4% (AVirgo,ALIG02,ALIG03) for NS-NS, from +3% 
(ALIGOl) to +12% (AVirgo) for BH-NS. 

Thus, it is clear that Approx 1+2+3, which has been 
used to estimate the detection rates of second generation 
detectors in p!8H2T] , can result in an overestimate or un- 
derestimate of the rates depending on the details of the 
population and on the interferometer which is considered. 



IV. CONCLUDING REMARKS 

In this paper we discuss the procedure to compute the 
detection rates for coalescing compact binaries, to be de- 
tected by the advanced detectors expected to be in op- 
eration in a few years. We quantify the effect of vari- 
ous approximations adopted in the literature to estimate 
these rates, and show that they introduce a bias which 
depends not only on the approximation which is used, 
but also on the statistical properties of the considered 
binary population (BH-BH, NS-NS, or BH-NS). These 
properties depend on the physical inputs introduced in 
the population synthesis codes used to generate the pop- 
ulation, and these may vary significantly from one model 
to another. An example is shown in Figure [2J where the 
difference between the mass distribution of BH-BH bi- 
naries for model A, which assumes solar metallicity, and 
model B, based on subsolar metallicity, is apparent. 

Thus, if we want to associate a detection rate to a 
given population of binaries, derived on the basis of some 
physical assumptions, we need to use the fully consis- 
tent procedure described in section [Hj indeed, the use of 
approximations as those described in section |III[ would 
not allow to disentangle the effects of the physical inputs 



from those of the assumed approximations. Finahy, we 
give an estimate of the detection rates for the Einstein 
Telescope, ET, using two population models obtained by 
the updated version of the St art rack code, which are 
publicly available. These rates have been obtained using 
the no approximation procedure, which is the only one 
applicable to ET, due to its high sensitivity. 
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